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Abstract 

Self-energy effects of Raman phonons in isotopically disordered semiconduc- 



^ ■ tors are deduced by perturbation theory and compared to experimental data. 

^D ' In contrast to the acoustic frequency region, higher-order terms contribute 

g: .^^.^.^^^.^^^.^...^^. 

ric dependence of the seff-energy of a binary isotope system mi_xMx on the 
concentration of the heavier isotope mass x can be explained by taking into 
account second- and third-order perturbation terms. For elemental semi- 
'^ \ conductors, the maximum of the self-energy occurs at concentrations with 



0.5 < X < 0.7, depending on the strength of the third-order term. Reasonable 
approximations are imposed that allow us to derive explicit expressions for 
r> \ the ratio of successive perturbation terms of the real and the imaginary part 

of the self-energy. This basic theoretical approach is compatible with Raman 
spectroscopic results on diamond and silicon, with calculations based on the 
coherent potential approximation, and with theoretical results obtained using 
ah initio electronic theory. The extension of the formalism to binary com- 
pounds, by taking into account the eigenvectors at the individual sublattices, 
is straightforward. In this manner, we interpret recent experimental results 
on the disorder-induced broadening of the TO (folded) modes of SiC with a 
^^C-enriched carbon sublatticcun 



Typeset using REVTgX 



I. INTRODUCTION 

The mass-fluctations (of a single element) induced by isotopic disorder can be described 
by the weighted n-th moments of the relative mass differences with respect to the average 
mass m = ^^ Cirrii of the virtual crystal defined as 

i ^ 

where i runs over the different isotopes with masses rrii and concentrations Cj. According 
to Eq. (P, f?! = within the virtual crystal approximation (VGA), i.e. if m is used as the 
atomic mass. 

For a two-isotope system mi_a.Mj. with M being the heavy mass and x = cm, Eq. (|1]) 
reads 

gn{x) = 1 - X ^^ + ^ ^^ • 2 

\ m J \ m J 

Equation (0) can be rewritten as follows, using Am = M — m 

gn{x) = {m/Am + x)~^ 

xs(l-a;)(s"-^ + (-l)"(l-a;)"-'). (3) 

For the elements under consideration here, the first factor of Eq. (^) leads to a weak x- 
dependence that slightly breaks the odd or even symmetry of gn{x) with respect to a; = 0.5. 
Neglecting this dependence (i.e. for x <^ m/Am) we can write 

(m/Am + x)^" ~ (Am/m)"'. (4) 

The lowest-order moments of Qn thus read (see Fig. |l]): 

g2 — (Am/m)^ x {1 — x) (5a) 

gs ~ (Am/m)^ x (l — x) (2x — 1) (5b) 

(74 ~ {Am/m^Y x (1 — a;) (1 — 3x(l — x)) (5c) 

g^ ~ {Am/mf x{l-x) {2x - 1) (1 - 2x(l - x)) (5d) 



Hence, the contribution to the self-energy is even (odd) with respect to x = 0.5 for even 
(odd) values of n. In particular, g2 represents a parabola with the maximum at a; = 0.5, 
while (73 is odd with respect to a; = 0.5, i.e. (73(0.5) = 0. 

The disorder-induced phonon self-energy is defined as 

ndis(^) = Adis(^) - 1 rdis(^)/2, (6) 

where Adis denotes the disorder-induced shift, and Fdis the disorder- induced broadening 
corresponding to the full width at half maximum (FWHM). The self-energy can be decom- 
posed into individual contributions that are proportional to the various moments of gn with 
n = 2,3,4,... 

ndis(t^) = A2x{l- x) 

+A^x{l-x){2x- 1) 

+AiX {I - x) {I - ?,x{l - x)) 

+ ... (7) 

In order to estimate the contributions of higher-order perturbation terms, Eq. ([Zp is fitted 
to experimental and theoretical data of group-IV semiconductors (see Sec. |T| for diamond. 



Si, Ge, a-Sn and Sec. |3 for SiC). In Eq. (|^) the fit parameters Ai = Aj a + ^A,r are 
complex numbers; in practice, we perform individual fits to Adis and Fdis, thus obtaining the 
two sets of fit parameters Aj^a and Aj^r- In the following section we show that, in general, 
the asymmetries with respect to x = 0.5 are not identical for Adis and Fdis. 

For a characterization of the contribution of the third-order perturbation term with 
respect to the second-order term, we define the ratio 



'~ A,,, 
with s = A, F. For the maximum self-energy Xmax,s we find 



(8) 



W. = ^±HJ3 + ;^-^|, (9) 



where the positive (negative) sign corresponds to r > (r < 0). 



II. PERTURBATION THEORY 

In this section, explicit expressions for the disorder-induced self-energies of a two-isotope 
elemental crystal are derived by perturbation theory. For the calculations we use Ildis = 
n = A + zr with A = A2 + A3 + . . . = Adis and r = Ta + Tg + . . . = -rdis/2. The 
zero-temperature phonon propagatoia D(uj, cUi) for an isotopically pure crystal is the sum of 
the propagators D^{u,u!i) and D_{u!,uJi) for positive and negative frequencies, respectively, 

D = D+{uj,uJi) + D_{uj,uJi) (10a) 

1 1 

+ (10b) 



UJ — UJi — I'-f —UJ — OOi — IJ 

= 2 ^2' , > (lOc) 

UJ^ — Uf — l2u!i'J 

where 7 — i> if anharmonicity is neglected. Ui denotes the absolute value of the phonon 
frequency. Since we are interested in the self-energy of the Raman mode, to which the 
largest contributions stem from close-lying optic phonon states, the propagator D^ can be 
neglected because |Re{D_}| ^ |Re{D+}| and |Im{Z}_}| ^ |Im{D+}|. In the following we 
thus set D = D^ with 

D+{uj,Ui) = — (11a) 

UJ — LUi — I'J 

+ ^7 X. , 9 (lib) 



{uj - uJiY -F 7^ {uj ~ LOiY + 72 
= Pi{u) + iU{u). (lie) 

Dj^iuo^uoi) consists of Pi{uj), with a pole at Ui, and a Lorentzian Li{uo). For the complex 
self-energy of the Raman phonon of a semiconductor with diamond structure, nth-order 
perturbation theory yielddil 

rr ( \ 3n{x) / 1 ^ "~^ 
lln(a;,a;) = uj ——- 



N oJiDi ■ ujjDj ■ ■ ■ uJnDn (12a) 



X 

i,j,...,n 



2" \ 6N, 



k?"'^*) 



'^hn^>."'-Di (12b) 



with Dn = D{uj,uJn) and Nc the number of unit cells. The second-order term {n = 2) 



contains the real part 



A2(a;,x) = ^^^ ( 7^J2^r .,^,J\. ] (13b) 



and the imaginary part 



r.(...) = *MJ 1 53^.^ ^_^ . (130) 

4 \ bi\c [^ — ^i) + 1 J 

For 7—^0, Eq. ( |13cD simplifies to a sum over 5-functions, which represents the one- 
phonon density of states (DOS) pi{uj), 

- u^, ^) = ^ ^' (^ E ^ ^(^ - ^o] (13d) 

= ^92{x)uj^piiuj). (13e) 

If A2(u;) and T2{uj) do not depend strongly on u near the Raman phonon at uq, then 
— r2(co'o) of Eq. ( p,3e| ) represents the half width at half maximum (HWHM) of the disorder- 
induced broadening. 

The third-order self-energy can be written as 

TT f \ 93{x) f I ^"^ 
li^iijj.x) = — - — UJ 



c 
X 

i,3 



^ uJiUjj{Pi + iLi){Pj + iLj), (14a) 



yielding the following real and imaginary parts 

A r ^ 33{x) f I ^^ 

A3[UJ,X) = — - — UJ 



X Y^ ^i^J (PiP] - ^i%) ' (14b) 

r3(a;,x) = ^^a;(^j ^UiUjPiLj. (14c) 



The corresponding fourth-order terms read 



A4(u;, x) = —— uj 



1 



16 " \6Nc 

i,j,k 



(15a) 



1 4(0;, x) = , „ UJ 



1 



16 \6N^ 

X 2_^^i^j^k{^PiPjLk — LiLjLk). 
i,j,k 



(15b) 



The ratios of the third-order to the second-order perturbation terms for the real and 
imaginary part of the self-energy is simply related to the disorder- induced shift A2 and 
broadening r2: 



A3 ^ gsjx) 1 
A2 2(72 (x) 6N, 



Yu^iPi 



J2 ^iL, 
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(16a) 



x=xo 



(16b) 



The ratios r^ (Eq. (|^)) for a certain phonon mode at ujq can be determined from the 
dilute limits x — > {0, 1} and the experimental shifts and broadenings for xq = 0.5. Under 
the approximation of Eq. (^), the expressions of Eq. ( |16a| ) and Eq. ( |16b| ) simplify to 

^ m 1 



rr — 16 



Am u; 
m 1 



A2 

[^2]^.„=0.5 



a;o=0.5 



(17a) 
(17b) 



Am UJ 

For phonons of crystals whose |A2(cc'o,Xo)| and |r2(ct;o,a;o)| are comparable, such as for 
the Raman mode of diamond or the TO Raman modes of SiC, the asymmetry of A{ujq, x) 
also depends on r2(co'o, Xq) and r^/r^ oc 2A2/(A2 — Eg) {rp/r^ ^ 8/3 for the Raman mode 
of diamond) . 

For phonons with |r2(ciJo)| <^ |A2(ciJo)|, such as the Raman modes of Si, Ge, and a-Sn, 
Eq. (|17a|) can be further simplified. In these cases, ta oc A2 and rr/rA — 2. The fact 



6 



that, generally, the contributions of higher-order terms to the self-energy are not equally 
distributed among the real and the imaginary part (r^ 7^ r-p) means, that the asymmetries of 
A{ujq, x) and T{ujo, x) have different shapes resulting in different concentrations of maximum 
self-energies (xmax.A 7^ a^maxj)- The ratios of the next successive-order terms read 

A^~AA^~2A^~2A^- ^ ^ 

Figure |^ shows the dependence on r^ of the concentration Xmax, at which the self-energy 
has a maximum, if |r2(ct;o)| ^ |A2(u;o)| (valid for Si, Ge, and a-Sn). 

General expressions for the ratios of successive-order perturbation terms for phonons 
with |r2(u;o)| ^ |A2(u;o)| can be easily found using 






= r/"(cos(n^) + i sin(n^)) ~ r/"(l + inO), (19) 



where rj = \D\ and 9 = Arg{D}. They are: 

^n+l _ 9n+l 



— V UiPi 1 



r„ 2gn n 

_ n An+i _ n A3 
n — 1 A„ n — 1 A2 



(20a) 

(20b) 
(20c) 



According to Eq. (|17a)) , negative values for r^ are possible if |r2| exceeds IA2I, but this does 
not occur for the Raman phonons of elemental semiconductors. One could also obtain rp < 
(Eq. ( |17b|) ) for phonons with q 7^ that exhibit a negative disorder-induced shift.Q The latter 
originates from the fact that the shift is related to the Kramers-Kronig transformation of 
the broadening, which displays a change of sign of Adis(^) near maxima of rdis(to'). In such 
a case the asymmetry would be flipped with respect to a; = 0.5, i.e. the maximum of the 
self-energy would occur at a concentration Xmax < 0.5. 



III. COMPARISON WITH EXPERIMENTAL AND THEORETICAL RESULTS 

In this section, we display a compilation of the disorder-induced self-energies for the 
Raman phonons of elemental crystals (diamond, Si, Ge, a-Sn) which have been obtained 
either by Raman spectroscopy or from theoretical calculations by several research groups 
during the past decade. Raman studies that address the variation of the self-energy with the 
isotopic composition have been conducted for diamond&ii^ and Si.li3 The coherent potential 
approximation (CPA) has been employed for diamondOa and Sic3, while ab initio electronic 
structure based calculations have been performed for diamond and Ge.ao All sets of data 
are analyzed with respect to the asymmetries of the self-energies and compared to the results 
we obtain by perturbation theory. 

Figure § displays the Raman frequencies of diamond versus the ^^C-concentration. The 
points (open symbols) represent experimental values. The dashed curve represents the 
approximately linear dependence expected in the VGA. The upward curvature of the ex- 
perimental data (with respect to the VGA line) clearly demonstrates the existence of an 
isotopic-disorder-induced self-energy as emphasized by the solid line, which is a fit with 
Eq. d^) for n = 2, 3. The prefactors A2 and ^3 contain all factors independent of x. The 
ratios r^ = ^43/^2 of the fitted values of these prefactors are listed in Table | under "Raman 
experiments". Note that two points, which are represented by dashed triangles, have not 
been included in the fit since they deviate unreasonably from the other data. 

It is difficult to see with the naked eye in Fig. ^ the asymmetric behavior versus x, 
which may arise from third-order perturbation terms. The asymmetry appears, however, 
rather clearly when the difference between the measured (or the calculated) behavior and 
the VGA line is plotted, as shown in Fig. ^. In this figure, the solid line also represents the fit 
to all experimental data (except the dotted triangles), the dot-dashed line represents GPA 
calculations while the dotted line is a fit to the asterisks which indicate points obtained in 
the ab initio calculations.aa All data in Fig. ^ show a similar asymmetric behavior, with 
a maximum of Adis(a;) at x ~ 0.6. This figure allows us to conclude that the real part 



of the self-energy due to isotopic disorder is well understood for diamond, including the 
superposition of second-order and third-order perturbation terms. 

Similar degree of understanding has been reached for rdis(x) as shown in Fig. |^. The 
x-position of the maxima of Ajis and Fdis determined from the experimental data agree 
with those obtained by perturbation theory (Eq. (^) and also with the CPA and ab initio 
calculations (see Table ||). 

Concerning the other elemental semiconductors, detailed experimental results with suf- 
ficient values of x to reach quantitative conclusions of the type found for diamond, are only 
available for Si.li3 These data for Si are shown in Fig. ^ and Fig. ^ (filled circles) together 
with the results of CPA calculations (filled squares). The latter show for Adis(x) a clear 
asymmetry with a maximum at Xmax,A ~ 0.56. The quality and the number of the ex- 
perimental points are not sufficient to conclude that an asymmetry exists but they cannot 
exclude it either. The measured absolute values of Adis(a;) almost (not quite) agree within 
error bars with the calculated ones. 

The corresponding experimental values of Fdis(x) (see Fig. |^) are about a factor of two 
lower than the calculated ones, although both show the asymmetric behavior (xmax,r ~ 0.62) 
predicted by Eq. ( |17bD . The reason for the discrepancy between the calculated and measured 



Fdis is to be sought in the mechanism responsible for it in Si.li3 Within the harmonic approx- 
imation Fdis = for Si, Ge, and a-Sn, because the Raman frequency is at the maximum of 
the spectrum and thus corresponds to zero density of one-phonon states. The rather small, 
but non-negligible, observed value of Fdis results from the DOS induced at the F-point by 
the anharmonic interactions responsible for the linewidth of the isotopically pure crystals. 
Thus, the widths observed for Si, as well as for Ge and a-Sn, correspond to fourth-order 
(twice disorder and twice anharmonicity) and higher-order terms. Under these conditions, 
we have no guarantee that the terms included in an anharmonic CPA calculation suffice to 
describe the experimental data, and agreement within a factor of two between theory and 
the small experimental values of Fdis is to be regarded as satisfactory at this point. Using 
this argument too, similar values of Adis are to be expected for these elemental crystals. 

9 



Although isotopic-disorder-induced effects have not been studied so extensively in Ge and 
a-Sn, there are some experimental results available for x = O.SEj'lj that can be compared 
with perturbation calculations. The latter, once contrasted with the measurements for 
X = 0.5, provide reliable predictions for future experiments. The corresponding data are 
displayed in Table |, where we compared them, in the case of Ge, with those obtained by ab 
initio calculations.t^M We should add that the ab initio values of Adis for Ge calculated by 
Vast and BaronitJEJ agree with the CPA calculation on Si and also represent the available 
experimental data rather well (^z^'*^ *"«**« (^ a; = o.5) = 1.14 cm~^ compared to 1.18 cm~^ 
and 1.06 cm^^ obtained experimentally for A2 in Si and Ge, respectively). The ab initio 
calculations of Fdis for Ge, however, give values which are an order of magnitude larger than 
the experimental ones, a fact that can be attributed to numerical errors in the calculated 
DOS near the van Hove singularity which is found at the Raman frequency. 

According to the values listed in Table |I| one can clearly distinguish between two different 
behaviors in the magnitudes of A2 and r2. Diamond exhibits a relatively large disorder- 
induced shift, A2, whereas the other elemental crystals have much smaller values. This 
largely reflects the factor of u in Eq. ( |13bD . The difference is even more striking for r2, 



which nearly vanishes for Si, Ge and a-Sn. This can be explained also by the different 
orders of perturbation theory coming into play. The overbending of the phonon dispersion 
of diamond introduces large (second-order) isotopic disorder effects which are not present in 
the case of Si, Ge, and a-Sn because of their vanishing DOS at the frequency of the Raman 
phonon. Such overbending of the phonon dispersion is also encountered in SiC and will be 



discussed in Sec. [V. 



We also observe a slight but systematic upwards shift of Xmax,r of the CPA values for 
Si, compared to those obtained by perturbation theory. The points of the CPA calculation, 
which includes higher-order perturbation terms, are fltted with Eq. (|^) for n = 2,3 that 
only contains the second-order and third-order terms. As shown in Fig. 0, the higher-order 
terms (such as, e.g., g^) emphasize even more the asymmetric appearance of the self-energy, 
provided that the prefactors Ai are all positive. A flt with Eq. (^ for n = 2, 3, 4 to the CPA 

10 



values for the disorder-induced broadening of Si can be found in Ref. [ O . 



IV. EXTENSION FOR BINARY COMPOUNDS 

The total isotopic-disorder-induced self-energy for a phonon of a binary compoundtll can 
be expressed as a sum of the self-energy contributions of the individual sublattices n 

K 

^ (^E^^^(^'^^)l^'^(^*)l') (21) 

where \e^\^ denotes the square of the eigenvector, averaged over all q's corresponding to a 
given frequency. The second-order self-energy can be described in a way similar to that used 
for monoatomic semiconductors 

K 
K 

x(3^i:-.|eV.)P (^_J).^^. ) (22b) 

For 7 — > 0, Eq. ( p2b|) simphfies to 

where plioj) corresponds to the one-phonon partial DOS projected on the sublattice k. 
Note in Eq. ( |22c| ) that for crystals with diamond structure Eq. ( |13e| ) is recovered by setting 

|e'^H|2 = l/2. 

11 



The general expression for r^ can be quite complicated for binary compounds when both 
sublattices contain isotopic disorder. However, if there is only one isotopically disordered 
sublattice, the sum over k is replaced by a single term. In this case, it is straightforward to 
derive expressions similar to Eq. ( |1 7a| ) and ( |17b| ). 



m 1 

rA - 



/\m ijj |e(u;)p 



A2 J ^=0.5 



(23a) 



rr ^ 16 ---r, [A2],=o 5 (23b) 

Am uj |e(u;)p 

where e(u;) is the eigenvector component at the disordered sublattice for the phonon 
under consideration (e.g. a Raman phonon). Thus, the isotopic effects of a disordered 
sublattice in a compound is different from that for the corresponding monoatomic crystal. 

We apply Eq. (|23a]) and Eq. (|23b| ) to the Raman spectroscopic results on a variety of 
'^'^'Si^^Ci-x^^Cx polytypes, recently reported by Rohmfeld et a/ [|^. We have performed a fit 
with Eq. (1^) for n = 2, 3 to the linewidth of the transverse optic (TO) modes of the 6H-SiC 
polytype measured in Ref. [ 0] for ^^C-concentrations ranging from x = 0.15 to x = 0.40. 
These experimental results are displayed together with our fits in Fig. |^. The fit performed 
in Ref. [ 0] was based on the asymmetric curve obtained by CPA calculations for diamond.tl 
This curve was first fitted to the data points of the TO (2/6) mode and further adjusted to 
the TO(0) and TO (6/6) modes by multiplication with the constant scaling factor uj'^pi{uj) 
assuming constant eigenvectors (see Eq. (|22cj ). Instead, we have considered each TO mode 
separately and performed fits with Eq. (|^) for n = 2,3 in the same manner as for the 
elemental semiconductors in Sec. |II[ We used, however, parameters appropriate to SiC, 
not to diamond. In this way, we rigorously conclude that the behavior of T^is versus x is 
asymmetric. This fact cannot be derived from the data in Fig. 3 of Ref. [ which were 
obtained only for 0.15 < a; < 0.40. The latter can be fitted equally with either a symmetric 
or an asymmetric curve. 

Using Eq. ( |22c| ) for a single disordered sublattice, together with the assumption of a 
constant eigenvector |eqi;o(q)| = 0.84,li3 we have calculated the imaginary part T2{uj) for 6H- 



SiC and 3C-SiC using the DOS of Hofmann et al. [ |T^]. The corresponding real part, ^2{^) 
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is obtained by the Kramers- Kronig relations. The frequency dependence of A2 and r2 for 
6H-SiC, as well as r2 for zincblende SiC (3C-SiC) is shown in Fig. ^ for the region of the TO 
modes. The three parallel lines cut the x-axis at the frequencies of the TO(0), TO(2/6) and 
TO (6/6) phonons, and the crossing with A2{QH) gives the (approximate) disorder- induced 
shift of the corresponding Raman peak. 

Raman experimental results extracted from Fig. ^ are compared in Table || with those 
found by the perturbation scheme of Eq. (|23a|) and (|23b| ). The fit to the TO(0) mode is 



symmetric, but the large scattering of data points still allows reasonable asymmetric fits with 
similar curvature than for the other two modes. The two numbers given for the ratio ta of 
the TO(0) mode represent maximum and minimum values and correspond to 6H-SiC and 
3C-SiC, respectively. For the cubic polytype the DOS vanishes at the Raman frequencies 
in the harmonic approximation. Therefore, we performed a convolution of the imaginary 
part with the anharmonic broadening Fanh = 1-4 cm~^ [ ^ in order to estimate the finite 



value. Note that the phonon dispersion of the hexagonal polytypes, calculated in Ref. [ |T9 
exhibits an overbending of the TO branches, which mainly occurs around the K-point of the 
Brillouin zone. The relatively large value of F2 observed for the TO(0) mode compared to 
those of Si, Ge and a-Sn (see Table |) provides evidence of a non-vanishing DOS caused by 
an overbending of the phonon dispersion. The experimental results from Ref. [ |^ support 
therefore the DOS calculated with the bond charge model by Hofmann et al. [ |19|] . 

The striking feature in Table ||, which is added to the general behavior of the Raman 
phonons of elemental crystals, is the negative values of r^ obtained for the TO (2/6) and 
TO (6/6) phonons. This change of sign implies a reversed asymmetry, i.e. the maximum is 
now shifted to lower concentrations of the heavier isotope (x < 0.5). Since in Ref. [ the 
phonon frequency was used for the determination of the ^^C-concentration, no experimental 
information can be exctracted about the disorder-induced shift. Within our interpretation, 
however, a negative ta has been already observed indirectly for a two-phonon combination in 
isotopically tailored diamond by cathodoluminescence.c^ The fact that some absolute values 
of rA,r are much larger than unity is caused by a strong difference in magnitude between A2 
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and r2 or by a vanishing A2 (see Eq. ( P3a| ) and Eq. ( |23b|) ). The latter occurs at frequencies 



near the maximum of the phonon DOS where A(ti;) changes sign. 

V. CONCLUSIONS 

General analytic expressions are given which allow the calculation of the real and imag- 
inary parts of the self-energy of phonons due to isotopic disorder in elemental crystals (e.g. 
diamond-type semiconductors) to various orders of perturbation theory. The results are 
used to analyze the asymmetric behavior of these self-energies versus the isotopic concen- 
tration in crystals containing two different isotopes. The asymmetry is found to originate 
primarily from a combination of second-order and third-order contributions to the phonon 
self-energies and can be determined from the second-order perturbation terms without the 
need of additional parameters. These results are extended to binary crystals (e.g. zincblende 
or wurtzite-type semiconductors) and used to interpret quantitatively recent experimental 
data for SiC versus the isotopic composition of the carbon component. The possible asym- 
metric behavior of the self-energy in this material, which contrary to the statement in Ref. [ 
H does not follow from the available data for 0.15 < x < 0.40 only, is clarified on the basis 
of our analytical expressions. 
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FIG. 1. Mass-fluctuations gn{x) for various n-th moments (Eqs. (g)) of an elemental, i.e. 
monoatomic, crystal (or each sublattice of a compound) containing a binary isotope system, x 
denotes the concentration of the heavier mass isotope. 

FIG. 2. Dependence of the value of Xmax at which the maxium of A and F is found, including 
second-order and third-order perturbation terms (solid line) and also fourth-order terms (dashed 
and dotted lines). 

FIG. 3. Raman shift of isotopically disordered diamond. The open symbols represent exper- 
imental values [ 0-0] The dashed line indicates the harmonic scaling of the phonon frequency 
within the VGA [uj oc m^^'^). The solid line corresponds to a fit with Eq. (j^) for n = 2,3 to all 
experimental data, added to the VGA-scaling. The two experimental data points at x = 0.1 and 
X = 1, represented as dotted symbols, were excluded from this fit. 

FIG. 4. Disorder-induced shift of the Raman phonon of diamond as a function of the 
^^C-concentration. The open symbols are Raman experimental data [ P-jTl]] whereas the aster- 
isks correspond to ab initio calculations [ 13,14|. The solid line is a fit with Eq. (^ for n = 2, 3 to 



all experimental data except for the two points at a; = 0.1 and x = 1 indicated as dotted symbols. 
The dotted and dot-dashed lines represent the fits to theoretical values obtained from ab initio and 
GPA calculations, respectively. 

FIG. 5. Disorder- induced broadening of the Raman phonon of diamond as a function of the 
^^C-concentration. The filled circles have been obtained from the Raman data of Refs. [ ^-10| by 



taking into account the corresponding instrumental resolutions and subtracting the anharmonic 
broadening, Fanh ~ 2 cm~^ (FWHM). The solid line is a fit with Eq. (^) for n = 2,3 to these 
points. The dotted and dot-dashed lines are the corresponding fits to the values obtained from ab 
initio [|1S 14] and CPA [ |^,^ calculations, respectively. 
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FIG. 6. Disorder-induced shift of the Raman phonon of Si as a function of the 
'^'^Si-concentration [ 12|. The sohd line is a fit with Eq. (0) for n = 2,3 to the experimental data. 



The dot-dashed line represents the corresponding fit to the values obtained from CPA calculations. 

FIG. 7. Disorder- induced broadening of the Raman phonon of Si as a function of the 
^'^Si-concentration [ |l2|. The solid line is a fit with Eq. ^ for n = 2,3 to the experimental data. 
The dot-dashed line represents the corresponding fit to the values obtained from CPA calculations. 

FIG. 8. Disorder-induced broadening of the TO Raman modes of the 6H-SiC polytype versus 
the ^'^C-concentration of the carbon sublattice. The data are taken from Ref. [|2[. The solid lines 
represent fits with Eq. (0) for n = 2,3 to the data points that correspond to the TO(0), TO(2/6), 
and TO (6/6) phonon modes. 

FIG. 9. Disorder-induced shift A2{uj) and broadening r2(u;) for 3C-SiC and 6H-SiC. We obtain 
r2(ci-') from the DOS that has been calculated by Hofmann et al. [|l^ using the bond-charge model. 
A2(w) is determined by a Kramers-Kronig transformation. 
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TABLES 
TABLE L Characterization of the asymmetry of the phonon self-energy with respect to a; = 0.5 

through the ratio r^ = ^3^3/^2,3 and the concentration Xmax.s with s = A,T. These quantities 

are displayed for each elemental crystal that is realized as a binary isotopic system rrii-xMx, with 

separate rows for the corresponding real and imaginary part. The Raman data are compared 

with theoretical results that are derived from perturbation theory (PTh, this work), the coherent 

potential approximation (CPA), and ab initio calculations. The values of the self-energies (A2,r2) 

are for x = 0.5. 



Raman experiments 



PTh 



CPA 



ab initio'^ 



A2 (cm ^) TA Xmax.A ^'A 2;max,A ^A Xmax,A TA Xmax.A 

mi_xMx ^ r---U 



-To (cm 



rr a^max.r rr Xmax.r ''r 2;max,r ''r 2;max,r 



12 






ga 


0.31 


0.57 


0.33 


0.58 


3"^ 


0.61 


0.62 


0.88 


0.66 



0.35 0.58 0.22 0.55 
0.66 0.63 0.53 0.61 







1.18" 


«0 




0.50 


0.26 


0.56 


0.25 


0.56 






28Sii_ 


x=^°Si. 


0.03'' 
1.06'= 


0.53 




0.61 


0.52 
0.33 


0.61 

0.58 


0.64 


0.63 


0.30 


0.57 


70Gei. 


-J'Ge^ 


0.03'= 
1.8'^ 








0.66 
0.67 


0.63 
0.63 






0.74 


0.64 


ii^Sm 


-x^^^Sn. 


0.02'^ 
0.7"^ 








^ 1 
0.41 


^0.67 
0.59 










ii^sm 


-xi24Sn. 


0.02'^ 








0.82 


0.65 










''Values taken as 


average from Refs. 


[|- 


-0. 














''Ref. 


0- 






















'^Ref. 


m- 






















^Ref. 


m- 






















^Refs. 


[13,14|. 
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TABLE II. Asymmetry of the phonon self-energies of isotopically disordered SiC, characterized 
by the ratio r^ = ^3^5/^2,3 and the concentration Xmax.s- Experimental results from Raman 
spectroscopy of the TO-phonons of 6H-SiC [ y] are compared with calculations using perturbation 
theory (PTh, this work). The values of the self-energies (A2,r2) are for x = 0.5. 



Raman experiment 



PTh* 



A2 (cm ^) TA 



:,A A2 (cm-1) 



natQ;12 



si^^ci_x^^a 



-r2 (cm 1) rr Xmax.r -r2 (cm ^) 



rr 



2^ in ax, A 
^max,r 



TO(0) 



TO(2/6) 



TO(6/6) 



iO.55 



i9.7 



0.71 



0.66 



iO.5 



0.64 



0.63 



1.7 


0.28-0.21 


0.57 


0.12-0.90 


0.58 


0.62 


2.2 


-0.75 


0.36 


3.8 


0.76 


0.64 


1.5 


-6.0 




7.3 


0.53 


0.61 



* Note that for negative values the asymmetry is flipped with respect to x = 0.5, i.e. the maximum 
of the self-energy occurs at x < 0.5. 
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